Astronomy & Astrophysics manuscript no. rstars 


© ESO 2009 


September 23, 2009 





The chemical composition of carbon stars. The R-type stars 

O. Zamora 1 , C. Abia 1 , B. Plez 2 , 1. Dominguez 1 , and S. Cristallo 13 

1 Departamento de Fisica Teorica y del Cosmos, Universidad de Granada, 18071 Granada, Spain 
■ e-mail: zamora@ugr.es 

q ■ 2 GRAAL, Universite Montpellier II, CNRS, 34095 Montpellier Cedex 5, France 

O ■ 

INAF, Osservatorio Astronomico di Collurania, 64100 Teramo, Italy 
C/3 ; ABSTRACT 



(N 



Aims. The aim of this work is to shed some light on the problem of the formation of carbon stars of R-type from a detailed study of 
their chemical composition. 

Methods. We use high-resolution and high signal-to-noise optical spectra of 23 R-type stars (both early- and late-types) selected from 
the Hipparcos catalogue. The chemical analysis is made using spectral synthesis in LTE and state-of-the-art carbon-rich spherical 
model atmospheres. We derive their CNO content (including the I2 C/' 3 C ratio), average metallicity, lithium, and light (Sr, Y, Zr) 
and heavy (Ba, La, Nd, Sm) i-element abundances. The observed properties of the stars (galactic distribution, kinematics, binarity, 
(-H ' photometry and luminosity) are also discussed. 

Results. Our analysis shows that late-R stars are carbon stars with identical chemical and observational characteristics than the normal 
(N-type) AGB carbon stars. In fact, the s-element abundance pattern derived can be reproduced by low-mass AGB nucleosynthesis 
models where the 13 C(a, n) ls O reaction is the main neutron donor. We confirm the results of the sole previous abundance analysis of 
early-R stars by Dominy (1984), namely: they are carbon stars with near solar metallicity showing enhanced nitrogen, low 12 C/ 13 C 
ratios and no ^-element enhancements. In addition, we have found that early-R stars have Li abundances larger than expected for 
post RGB tip giants. We also find that a significant number (~ 40%) of the early-R stars in our sample are wrongly classified, being 
probably classical CH stars and normal K giants. 

Conclusions. On the basis of the chemical analysis, we confirm the previous suggestion that late-R stars are just misclassified N-type 
carbon stars in the AGB phase of evolution. Their photometric, kinematic, variability and luminosity properties are also compatible 
with this. In consequence, we suggest that the number of true R stars is considerably lower than previously believed. This alleviates 
the problem of considering R stars as a frequent stage in the evolution of low-mass stars. We briefly discuss the different scenarios 
proposed for the formation of early-R stars. The mixing of carbon during an anomalous He-flash is favoured, although no physical 
mechanism able to trigger that mixing has been found yet. The origin of these stars still remains a mystery. 
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C> 1 . Introduction chemical peculiarities can be expl ained by the sec ond scenario 



as a consequence of mass transfer dHan et all ! 1995). 



Carbon stars are easily recognisable by the presence of absorp- The R stars seem to be very common among t he giant carbo n 

53 ' tion bands of C 2 and CN at near visual wavelengths and are stars, accounting 10 times more than the N stars (tBlancol[T965l). 

chemically characterised by C/O > 1 in their envelope. The Acco rding to the general catalogue of carbon stars dStephensonl 

Henry Draper classification divided carbon stars into the N and I l973l). R-type stars m ay amount to ~ 30 % of all the carbon stars. 

R groups on the basis of their spectral features. The N or nor- Bergeat et all d2002al) suggested that this fraction might be even 

mal carbon stars, show a very strong depression in the spectrum larger. This figure is very important since they may represent a 

at A < 4500 A, while the R stars seem to be warmer and the stage of evolution that is available to an appreciable fraction of 

blue/vi olet reg i on of the spectrum is usually accessible to obser- stars, and are not the result of anomalous initial conditions or sta- 

vations. lShand d 1 9281) split the R class into R0 to R8, where RO-4 tistically unlikely events. Further, their velocity dispersion and 

(hot-early-R stars) are warmer and equivalent to the K-type gi- position in the Galaxy indicates that early-R stars are members 

ants, and R5-8 (cool-late-R stars) are equivalent to M stars. Since of the galactic thick disk while late-R show kinematic p roperties 

C/O < 1 almost everywhere in the Universe, the carbon excess rather similar to the thin disk stellar p opulation (e.g. Sanford 

must be a result of stellar nucleosynthesis within the star itself ll944tlEggenlll972tlBergeat et al.ll2002bl) . The luminosities of R 

or in a binary companion. The N stars are understood in the for- stars (at least the early- R) are know n to be too low to be shell he- 

mer scenario: they are luminous, cool stars with alternate H- and Hum burning star s (e.g.lScaloll976l). In this sense, a fundamental 

He-burning shells, that owe their carbon enhancement to the mi- step was made by Knapp ^t alJd200 l[) wh ose re-processing of the 

xing triggered by the third dredge-up during the asymptotic giant Hipparcos data dPerrvman & ESAll 19971) . located R stars in the 

bran ch (AGB) phase in the evolution of low-mass (< 3 M Q ) stars H-R diagram at the same place than the red clump giants and 

(e.g. llben & Renzinllll983l) . Other carbon stars like the classical concluded that they were He-core burning stars or post helium 

CH-type are observed to be members of binary systems and their core stars. Another fundamental property of these stars is that no 
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R star has been found so far in a binary system (McC lurell 19971) . 
a statistically unlikely result. 

The only chemical analysis by iDomin showed that 

early-R stars have near solar metallicity, low (< 10) 12 C/ 13 C ra- 
tios, moderate nitrogen excess and no ^-element enhancements. 
No chemical analysis of late-R stars exists up to date. This find- 
ing for early-R stars is also in sharp c ontrast with that found in 
N-type carbon stars ( I Abia et al]|2002l) . So, the problem is how a 
single giant star not luminous enough to be on the AGB phase 
can present a C/O > 1 at the surface. The favoured hypothesis 
so far is that the carbon produced during the He-flash is mixed 
in some way to the surface. However, standard one-dimensional 
He-flash models do not predict the mixi ng of carbon-rich mate- 
rial fr om the core to the stellar envelope dHarm & Schwarzschildl 
1966). A lively discussion about this subject was held follow- 
ing die_^re^v]£gestionof mixin g at the He-flash proposed 
by ISchwarzschild & Harm! d 1 962b. Mod els in which rotation 
was parametrised dMengel & Gross|[l976l) lead to off-center He- 
ignition but not to conditions favouring mixing. On the other 
hand, models in which the location of the He-ignition was ad hoc 
move d to the outer part of the He-core dPaczvnski & Tremaind 
|1977|) produced the desired mixing. The general conclusion be- 
ing dDespainll 19821) that hydrostatic models do not produce mix- 
ing at the He-flash. The situation is differe nt for Z=0 models 
dFuiimoto et al.lll990l:lHollowell et al.lll990h but in this case the 
mixing is directly linked to the lack of CNO elements in the 
H-shell and it can not be extrapolated to more metal-rich mod- 
el£]. Recently, two- and three-dimensional hydrodynamic sim- 
ulations of the core He-flash in Population I stars h a ve been 
perfo rmed dLattanzio et all 120061: iMocak et ai1l2008allb1: iMocakl 
2009). None of these simulations lead to substantial differences 
with respect to the hydrostatic case. Thus, no physical mecha- 
nism able to trigger the required mixing of carbon during the He- 
flash in single stars has been found up to date. Binary star merg- 
ers have been invoked to induce such mixing. This, b y passing, 
might explain why no R star is found to be binary. Ilzzard etafl 
d2007l) investigate statistically possible channels for early-R star 
formation by a binary merger process. They found many pos- 
sible evolutionary channels, the most common of which is the 
merger of a He white dwarf with a hydrogen-burning red giant 
branch star during a common envelope phase. However, it is far 
from clear if such a merger might lead to the mi xing of carbon 
to the surface during the He-flash (Zamor a 200.3). 

In order to shed light on the problem of the origin of these 
stars, we present here a detailed chemical analysis of late- and 
early-R stars using high quality spectra. This is a further step in 
the study of the chemical composition of fiel d giant carbon stars 
in the Galaxy. In Paper I (Abia & Isern 2000) we studied J-type 
stars and in Papers II and III dAbia et al.ll200lLl2002T) the N-type 
stars. In the next section we describe the main characteristics 
of the sample stars and the observations. Sect. [3] describes the 
determination of the stellar parameters and chemical analysis, 
and Sect.|4]discusses the results and draws some conclusions. 

2. The data 

2.1. Observations 

We selected 23 galacti c R-type stars fro m the sample of car- 
bon stars compiled by Knap p et al.l d2001h with measured para- 

1 We remind that from the analysis by Dominy ( 1984), R stars have 
[Fe/H]~ 0.0. In this paper we adopt the usual notation [X/H]=log 
N(X)/N(H)*-log N(X)/N(H) Q , where log N(H)= 12 is the hydrogen 
abundance by number. 



flaxes according to Hipparcos dPerrym an & ESA 1997). The 
stars were observed with the 2.2 m telescope at CAHA ob- 
servatory (Spain) during March 2003 and July/August 2004. 
Spectra were obtain ed using the FOCES echelle spectrograph 
(Pfei ffer et al.l ll998). The spectra cover the wavelength region 
from AA ~ 4000 - 10700 A in 90 orders with full spectral cov- 
erage, a resolving power of R~ 40 000 in most of the spectra ex- 
cept those obtained in the last period of observation (R~ 20 000), 
because of the 2x2 CCD binning made to increase the signal-to- 
noise ratio in the spectrum of the fainter stars. The typical S/N 
ratio reached in the spectra was larger than 100 in the red orders 
and lower than this value below A ~ 5000 A. In fact, in some 
stars the blue orders were not useful for the abundance analysis. 

The spectra were reduced using the echelle task within the 
IRAF package following the standard procedure. When several 
images of the same object were obtained, they were reduced 
independently and finally added using the IRAF task scorn- 
bine. The wavelength calibration was always better than 0.03 
A. Finally, the individual spectral orders were normalised to a 
pseudo-continuum by joining the maximum flux points and, in 
some cases, corrected with the help of the synthetic spectra (see 
details of the method in Paper I). We estimate an uncertainty 
in the continuum location of less than a 5%, although it may 
be larger for the spectral orders severely affected by molecular 
bands, particularly in late-R stars. 

2.2. Spectral classification 

One of the most serious problems when dealing with carbon 
stars in general is their ascription to one or another spectral type. 
This is because the spectral classification is done frequently by 
using low resolution spectra, which does not permit the detec- 
tion o f relevant atomic lines due to blends w i th molecular bands 
(e.g.lCannon & Pickerindl918t]ShanJl928l:lKeenan & Morganl 
1941; lVandervordl958HKeenanll993tlBarnbaum et al.ll 19961) . In 
this sense, studies on early-R stars are frequently contaminated 
with CH -type stars, whose spectral d istribution is very similar 
(see e.g. IWallerstein & Knapplll998l) . In fact, the most effec- 
tive way of distinguishing between these two spectral types is 
to compare the intensity of some metallic lines (less prominent 
in CH-type stars) and, mainly, the intensity of s-elem ent lines 
which are well known to be enhanced i n CH-type stars dKeenanl 
119931: lAbia et al.ll2003t lGoswainill2005l) . An additional problem 
concerns the derivation of the temperature subtype. In general 
the intensity of C-bearing molecular bands in carbon stars de- 
pends not only on the effective temperature but also on the ac- 
tual C/O ratio in the atmosphere. In the case of early-R stars this 
introduces an extra uncertainty in the assignation of the tempera- 
ture subtype compared with normal (O-rich) G- and K-type stars 
of similar effective temperature. For the late-R stars, there is no 
tight correlation between the spectral distribution and effective 
temperature due to the presence of very intens e molecular ab- 
sorptions, similar to N-type stars dKeenanl fl993). TableQ]shows 
the spectral classification of our stars from different sources in 
the literature. It is evident that many stars have been classified in 
very different ways because of the problems cited above. 

This confusion in the spectral classification of many R-type 
stars has motivated the us e of others criteria such as the pho- 
tometric one. For instance, Knap p~et~ai1 (1200 ll) use the (V - K) 
colour index to distinguish between early and late-type R stars 
according whether (V - K)o < 4 or (V - K)o > 4, respectively. 
Since many colour indexes might be affected by the variability of 
the star (at least for the late-R) and also by the presence of strong 
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Name 




Var. Type Date 


Exp. 


Spectral Type 


HIP 


BD 


HD 


GCVS 




(*) 


(s) 




81 


-03°1873 


57 884 


V / JO 1V1U11 


Lb 


1 


^000 
jUUu 


N (Q\ R8 n 9 s ! 




+24° 1686 


59643 


LNV^ vJclil 


SR 


1 




Ky yizf, ko \yf, Ko V^-^/ 


7Q 1 1 8 
jy 1 lo 


- 


- 




- 


1 


J)OWW 


R9 ("1 


ztzt81 9 


- 


78 278 




- 


1 




DA /Q\ DC n 


C7 Q7 

J j ojz 


+41°2150 


- 




- 


1 


J)OWW 


KU ^UJ, V^rl-llKe (10J, V^rl \Z.) 


Jo / oO 


+71°600 


- 




- 


1 


JOWW 


KZ 1.1 V^rl-llKe ^1 / J, V^rl yZ J 


A 9 Ztf)1 
OZ ^fUl 


+04°2651a 


Ill 166 


PT T Vir 

ku vir 


Mira 


1 


JOWW 


p; (Q\ p; 

Kj l^y ) 7 Kj 1.1 -V 


A 9 QAA 


- 


112127 




- 


1 


Z^WW 


KU 1.1.5J, Kl \L^r) 7 KZ (il, o, 1UJ, I\J ^ / ) 7 Kj (,1 J 


to ARQ 
Oy KJoy 


- 


123 821 




- 


1 


JOWW 


no ct\ r,a n 1"> R9 czn R9 (i\ 

Oy ( / ), Uo ( 1 1 J, KZ (4 J, KZ ^ 1 J 


7zt 89A 
1 1+ ozo 


+30°2637 


- 




- 


1 


9700 
Z / WW 


RO H 9 s ! 
RU 1,1 Z ^ 


89 1 8zt 
oZ 1 o^f 


+23°2998 


- 




- 


3 


jOWW 


RO H ^"1 R9 ff\ 
KU 1,1 J,), KZ 1^1 J 


8/1 9AA 
o^f ZOO 


+42°2811 


156074 




- 


2 


Z^fWW 


RfWI 1 M\ R1 R9f1"\ 
KU 1, 1 1 , 1 / ), Kl ^JJ, KZ ^ 1 J 


OJ / J\J 


+02°3336 


- 




- 


3 


J>OWW 


R9 n CIA lilrn ('I 7"\ (A\ Mzt f 1 ~\ 
KZ ^10 J, V^rl-llKe 1,1 / IN J \ ' /, ^ ' 


oo yz / 


+17°3325 






- 


2 


j www 


RO ("1 9~\ 
KU ^ 1 Z J 


87 Afn 




1 A 7 878 
10j OJO 






J 


J'+WW 


p; (Q\ p« /i c\ 
Kj {y ) 7 KJ (^IJJ 


88 887 


+09°3576 


166097 


_ 




2 


3600 


R4(9),R5 (15) 


91929 


-13°5083 


173138 


RVSct 


Lb 


3 


3600 


R3 (9), R3 (15), R4 (1) 


94049 


-17°5492 


178316 






3 


5400 


R4 (12), R2 (9) 


95422 










3 


5400 


R5 (15) 


98 223 


-00°3883 


188 934 






3 


7200 


R8(12),R4(9) 


108 205 


+49°3673 


208 512 


LWCyg 


Lb 


3 


4200 


R3 (9),R2 (15), R8 (3) 


109158 






CT Lac 


SRa 


3 


5400 


R (6), N8 (3) 


113 150 




216649 






3 


3300 


R3 (9),R5 (15) 



* (1) 2003 March 12-13; (2) 2003 August 9-10; (3) 2004 July 1-4. 

R eferences for th e spectral ty pes: 

1: Barnbaum et al. 1 1996); 2: Bartkevicius 1 1996); 3: Ealitis et al. 1 2003); 4: Keenan 1 1993); 5: Keenan & Morgan 1 1941); 6: Lee & Bartlett i 1944); 7: McClure i 1970); 8: Morgan & Keenan 
119731); 9:|SanforJ 119441) ; 10: ISchild 119731) ; ll:ISchmitt 1 1971); 12: Shane 1 1928); 13: Stock & Welhau i 1959): 14: IUpgred 119621) : 15: IVandervonl 119581) : 16lYamashita (1973); 17: 
Yamashita i 1975). 

The following stars show evidence of binarity: 

HIP 36623 ( Carquillat & Prieur 2008), HIP 109 158 iMakarov & Kaplan 2003), HIP 53 832 iPlatais et al. 2003) and HIP 85 750 i McClu rel 19971) . 



molecular absorptions, we simply adopt here the criterion based 
on the effective temperature estimated in the chemical analysis 
(see Sect.O. We will see, that an effective temperature threshold 
of ~ 3600 K seems to be a good criterion to distinguish between 
early- and late-R star s. This will lead to a spectral classification 
in agreement with the Knapp ^taT] d200ll) 's criterion above, and 
with other characteristics that differentiate the two subtypes of R 
stars. 



2.3. Distribution in the Galaxy and kinematics 

That the R s t ars are galactic disk o bjects was recognised by 
lEggenl dl972l) . iBergeat et all d2002bl) calculated the space den- 
sity in the galactic plane of early-R stars and found that it is a fac- 
tor ~ 16 lower than for N stars; in fact, on average they are three 
times further from the galactic plane. Despite the small sample 
studied here, we reach the same conclusion: our early-R stars are 
located at large galactic latitudes, |b| > 30°. The distribution of 
the late-R stars is however, alm ost identical to that of the N-type 
stars (e.g. lClaussen et aflll987l) : they are very close to the galac- 
tic plane. These difference s in the galactic distribution of R stars 
were already known (e.g. ISanfordlll944t llshidal 1 1 960t [Rybski 
1 19721: IStephensonlll973l: iBarbaro & Dallaportalll974l) and are 
indicative of early-R stars belonging to the galactic thick disk 
while late-R stars belong to the thin disk. This of course implies 
a range for the masses and ages of the R stars: early-R stars must 
be of lower masses (< 1 M ) and older than late-R stars. 



Previous kin ematics analys e s perfor med by e.g. lDeanldl972l) 
and recently by Berge atet al.l (|2002b) show that the velocity 
dispersion of R-stars are typically larger by a factor of ~ 2 
than for N-type stars. IBergeat et al . (2002b) in particular, obtain 
<x = 42 - 54 km s~' in the direction of the galactic north pole 
for the stars belonging to their hot carbon grourfl where most 
of our early-R are included, and 23 km s~' typically for the cool 
variable group, which includes the late-type stars in our sample. 
This reinforce the conclusion that the differences between early- 
and late-type R stars are representative of two different stellar 
populations, as we already noted. According to t his numbers, 
and a pplying a standard age-velocity relation (e.g. lWielen et al.l 
1992), an age of ~ 3 Gyr and > 10 Gyr is obtained for late- and 
early-R stars, respectively. 

2.4. Binarity 

Table Q] indicates the stars with evidence of binarity obtained 
from radial velocity variations or other methods. These stars are 
HIP 53 832 and HIP 85 750, classified as early-R, and the late- 
R stars HIP 36623 and HIP 109 158. It is important to note 
that these two early-R stars have been also classified as CH- 



2 Bergeat et al. (2002b) define 14 photometric groups in order to clas- 
sify the carbon stars observed by Hipparcos into homogeneous classes. 
They define the hot carbon (HC) group which includes mostly early- 
R and CH-type stars, whereas the cool variables (CV) group includes 
N-type stars and a few late-R stars. 
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Table 2. Photometric data and luminosities 



Star 


V 


K 


J-H 


H-K 


Afro 


M hoX 


Ref. 


Late-R stars 


HIP 35 810 


9.01 


3.67 


0.83 


0.42 


-6.84 


-4.05 


1 


HIP 36623 


8.02 


2.95 


0.72 


0.38 


-7.23 


-4.18 


5 


HIP 62401 


11.98 


1.81 


1.85 


1.30 


-7.48 


-5.23 


7 


HIP 91 929 


9.75 


3.41 


1.30 


0.51 


-7.25 


-4.96 


1 


HIP 108 205 


9.23 


1.71 


1.36 


0.80 


-8.57 


-5.45 


4, 1 


HIP 109 158 


10.12 


2.57 


1.21 


0.75 


-8.90 


-5.65 


5, 1 


Early-R stars 


HIP 39118* 


10.41 


7.23 


0.79 


0.25 


-2.18 


0.47 


1 


HIP 44 812 


10.61 


6.50 


0.78 


0.35 


-4.73 


-2.02 


1 


HIP 53 832 


10.11 


7.60 


0.52 


0.13 


-2.83 


-0.81 


1 


HIP 58 786 


10.27 


7.51 


0.52 


0.15 


-3.06 


-1.02 


1 


HIP 62944* 


6.91 


4.16 


0.55 


0.10 


-1.27 


0.86 


3 


HIP 69089* 


8.68 


6.40 


0.46 


0.12 


-3.09 


-1.24 


1 


HIP 74826 


9.78 


7.39 


0.48 


0.12 


-3.31 


-1.38 


1 


HIP 82184 


9.10 


6.46 


0.57 


0.12 


-3.67 


-1.04 


1 


HIP 84266 


7.60 


5.11 


0.30 


0.23 


-3.09 


-0.87 


5 


HIP 85 750 


9.42 


5.15 


0.82 


0.24 


-5.85 


-3.25 


2, 1 


HIP 86927 


8.71 


6.14 


0.53 


0.14 


-1.69 


0.22 


2 


HIP 87603 


10.72 


7.92 


0.55 


0.18 


-2.24 


-0.02 


1 


HIP 88 887 


9.78 


5.13 


0.84 


0.36 


-5.90 


-3.14 


1 


HIP 94049 


10.39 


7.39 


0.55 


0.19 


-3.82 


-1.97 


1 


HIP 95422 


11.01 


6.74 


0.75 


0.26 


-4.81 


-2.57 


1 


HIP 98 223 


9.36 


5.15 


0.78 


0.31 


-6.27 


-3.78 


6, 1 


HIP 113150 


10.75 


7.97 


0.55 


0.22 


-3.45 


-1.26 


1 



Non dereddened VJHK magnitudes. 

References of photo metry : 

1: 2MASS ICutrietalJ (20p); 2: IDominv et alj UMS 
Neusebauer & Leishton 1 1969); 5: Noguc hieHtj] 119811) : 
fTJil); 7: Whitelock et al. 1 200(il). 



3: lElmsl (vWti: 4: 
6: iMendoza & Johnson! 



type stars, which are known to be all binaries. We will reach 
the same conclusion on the basis of their chemical composition 
(see Sect. 4.1). HIP 85 750 is a binary s ystem whose orbital pa- 
rameters were determined by McClure] d 1997b . The late-R star 
HIP 36 623 is a symbiotic star detected by time-variations of the 
ultra-violet continu um due probably to the presence of a white 
dwarf companion dJohnson et al.lH988t iBel czvriski et al. 2000; 
iMunari & Zwitterl |2002|) . Recently, ICarquillat & Prieurl d2008 ) 
estimate that the components of the system have masses of 2.5 
M Q and 0.6 M©. On the other hand, HIP 109 158 appears in 
the catalogue of Hipparcos astrometric bina ries with accelerated 
proper motions (Makarov & Kaplan 2005). For the remaining 
sample, there is no evid ence of binarity (in agreement with the 
study by iMcClurell 1997b or available study. 
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Fig. 1. Near infrared colour-colour diagram for the stars in this 
study. Triangles are early-R stars, fi lled circles l ate-R stars and 
empty circles the N-type stars from lAbia et al.l (1200 ll) . The re- 
gions typically pop ulated by CH- and N-type stars are indicated 
dTotten et alfeOOOh . 



and the relations by ICardelli et all dl989b . For the stars HIP 
39 1 18, HIP 62 944 and HIP 69 089, not included in the Bergeat 
et al. (20 02a) sample, we use the galactic interstellar extinction 
model bv lArenou et al.l (Q992) and the parallaxes by Knapp et al. 
(2001). In Fig.[TJthe position of our sample stars is plotted in the 
near infrared colour-colour diagram. We have marked with lines 
the re gions usually occu pied by CH- and N-type stars accord- 
ing to lTotten et al.l 12000). Some galactic carbon stars of N-type 
(Paper II) are included for comparison. It is clear from this fi- 
gure that early-R stars occupy the same region as CH-type stars, 
while late-R stars are redder and are located mostly in the region 
of N-type carbon stars. 

Another interesting photometric index is the K- [ 1 2 colour, 
which is considered a t racer of the dust surrounding the star 
(Jonssen &Knapdri 998). The formation of circumstellar dust is 
considered a measure of the stellar mass-loss and in consequence 
might give an indication of the evolutionary status of a given star: 
stars with no mass-loss or very low rates have K - [12] < 0.7. 
Our early-R stars fulfil this criterion and do not present sig- 
nificant mass-loss. At contrary, most of our late-R stars show 
K — [12] > 0.7 and would correspond to stars with moderate 
mass-loss rates, 5x1 0~ 8 to 5 x 1 -7 Mq y~ 1 , the extreme case be- 
ing HIP 62401 (K - [12] = 4.1) having a considerably larger 
rate, 10~ 6 Mq y -1 (the star in the upper right corner in Fig. [T). 
These number s are similar to those estimated in normal AGB 
stars (see e.g. Bussoetal . 20071). 



2.5. Photometry 

Table El shows the available VJHK photometry taken from 
Knap p et"ai] d200ll) . the Hipparcos catalogue or from the 
2MASS (Two M icron All Sky Survey) on-line data release 
dCutri et al.ll2003l) . Differences in the photometry found in the 
literature for a given star are typically below ~ 0.2 mag for the 
majority of early-R stars, but it can be higher for late-R stars 
due to their known photometric variability (see Table [T). All 
the photometric magnitudes were correct ed from interstellar ex- 
tinction using the A j values derived by Berge atet alj ((2002a) 



2.6. Luminosity 

As we mentioned above, the selected stars have all mea- 
sured trigonomet r ic p arallaxes according to Hipparcos 
( Perrvma n & ESAl 1 1997b . However, Hipparco s parallaxes 
for gi ant stars resulted to be not very accurate. Knap p et all 
d200ll) reprocessed the original Hipparcos parallaxes by im- 
posing the co ndition that they remained always greater than 
zero (see also iPourbaix & Jorissenll2000l) . In the sample of R 



3 [12] 
at 12 fim. 



-2.5 log (F12/28.3), where F12 is the IRAS flux density 
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stars in iKnapp et al.l d2001l) only 18 % of the stars have n/e(n) 
> 2, where n d enotes the tr i gono metric parallax and e{n) the 
associated error. lKnapp et al.l d200ll) perf ormed a re-reduction of 
the In termediate Astrometric Data (IAD. Ivan Leeuwen & Evansl 
1998) taking as parameter, in the x 1 minimisation, the logarithm 
of the distance instead of the distance itself. In this way, 
the derived parallaxes are always greater than zero. The true 
parallaxe s, i.e. the paral l axes free of all the possible biases 
(see e.g. lArenou & Luril [[999), 

were then estimated with a 
Monte-Carlo simulation by adopting diffe rent distributions for 
the ab scissa residuals of the IAD (see also Pou rbaix & Jorissenl 
2000). Following this procedure, these authors derived a true 
average absolute ^-magnitude for the early-R stars: M K „ - 
-2 + 1 , similar to the value -1.61+ 0.03 found bv lAlvesI (2000) 
for red clump giant stars near the Sun observed by Hipparcos. 
Another major study on the parallaxes of R stars was made by 
Bergeat et al. (2002a). These authors use the original Hipparcos 
parallaxes and correct them from all the observational biases 
(that may introduce an uncertainty u p to ~ 0.4 mag in the 
estimation of the luminosity, see e.g. Bergeat et al. 2002a| for 
further details), obtaining significantly different results. In fact, 
they find an average ^-absolute magnitude of Mk - -3.0 
for early-R stars. Correspondingly, they derive also bolometric 
magnitudes for R type s t ars ~ 1-2 mag brighter than those 
obtai ned bv IKnapp et al.l (1200 ll) . We believe that the parallaxes 
by Bergeat et al. (2002a) are more accurate since the parallaxes 
bv lKnapp et al.l d200ll) might undergo a bias when K > 8, whose 
effect is to overestimate the Mk (see their Fig. 2). On the other 
hand, a new analysis of the Hippar cos data has been performed 
recently by Ivan Leeuwenl (120071) . The absolute magnitudes 
computed with these new data for our stars (13 o f them) are in 
general in good agreement with those derived by Bergea fet al.l 
d2002al) . with a maximum difference below 0.5 magnitudes. 
Therefore, w e have adopted here the parallaxes derived by 
Bergeat et al. (2002a) to estimate the luminosity of our stars. 
For HIP 39 118, HIP 62 944 and HIP 69 089 (not included in the 
Berg eat et al.l 2002a sample, as previously noted) we use instead 
Knap p et al.l d200ll) parallaxes. 



Table [2] shows the absolute magnitude in the K band cor- 
rected for interstellar extinction (M^ ). The corr esponding ab- 
solute bolometric magnitudes (Mt, i) according to lBergeat et alj 
d2002al) are also shown. For the stars HIP 39 118, HIP 62944 
and HIP 69 089 we u sed the bolometric correction in K by 
ICosta & Frogell d!996l) . We note that the error in the derivation 
of M Ko and M\, Q i can be as large as ~ 1-1.5 mag, and it is 
dominated by the uncertainty in the trigonometric parallaxes. 
Nevertheless, we have constructed the H-R diagram of the stars 
in the sample (see next section for details in the derivation of the 
Jeff)- From Fig.|2]it is clear that the luminosities of late-R stars 
are close to the values expected for low-mass stars (< 2 M ) in 
the AGB phase. In Fig. [2]is also drawn the minimum luminos- 
ity (dashed line at M\, \ * -4.5) for a 1.4 M G st ar to become a 
therm ally-pulsing AGB carbon star according to iStraniero et alj 
(2003). This luminosity limit is quite model dependent (metal- 
licity, mass loss history, third dredge -up efficiency, treatme nt of 
the opacity in the envelope etc, see IStraniero et al] [2003). but 
roughly indicates that late-R stars have luminosities compatible 
with the AGB values. At contrary, early-R stars are clearly below 
the expected luminosity in this phase. 



-6 



-4 



TP-AGB, C-stor 


• 


t 




A 

▲ 
AA 


• 




- 

















3.70 3.60 3.50 3.40 3.30 

log T eff 



Fig. 2. H-R diagram for the sample stars. Solid triangles are 
early-R stars and circles late-R stars. Open triangles correspond 
to stars classified as early-R stars in the Hipparcos catalogue 
which have been reclassified here as K- or CH-type giants (see 
Sect. 4). 



3. Atmospheric parameters and abundance 
calculations 

We use the latest generation of MARC S spherically-symm etric 
model atmospheres for carbon stars dGustafsson et alJl2008l) . In 
this new grid however, there is no a complete set of carbon en- 
hanced models for T e g > 3500 K. For stars with effective tem- 
perature larger than this value, we use MARCS C-rich models 
computed by one of us (B. Plez). 

The estimation of the atmospheric parameters was made in 
a similar way to Papers II and III, and will be only briefly de- 
scribed here. We refer to these previous works for more details. 
We followed an iterative method that compares computed and 
observed spectra modifying T e g, gravity, average metallicity and 
the C/O ratio from given initial values until a reasonable fit to 
the observed spectrum is found. For the effective temperature, 
initial values were derived as the average temperature obtained 
from the calibr ations of (V - K)p , (H - K)o and (J - K)q vs. T e s 
according to Berge atet al.l d200ll) . For a given object, the typical 
dispersion in the T e s derived from the three indexes was +200 
K for the late-R stars and +500 K for early-R stars. However, in 
the case of the early-R stars, the theoretical spectra are so sen- 
sitive to a small variation of the T e s, that the uncertainty in the 
temperature is certainly much lower than ~ + 500 K. The spec- 
troscopic method of deriving the effective temperature is not rec- 
ommended in our case because of the blending with molecular 
lines, even in the warmer early-R stars. Further, the range in ex- 
citation energies of the apparently clean atomic lines available 
in our spectra is too narrow. Nevertheless, we checked that for 
the warmest early-R stars where the molecular contribution can 
be considered weak, we did not find any correlation between the 
estimated abundance of iron and the corresponding excitation 
energy of a number of Fe I lines (for instance in HIP 84 266, 
HIP 69 089 and HIP 39118). 

Gravity was estimated using the relation between lumino- 
sity, effective temperature and stellar mass. We adopted a mass 
of 1 Mq for all the stars in the sample. Although late-R stars 
are probably more massive (see Sect. 2.3), an uncertainty of a 
factor of two in the stellar mass translates to an uncertainty of 
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Fig. 3. Observed (filled circles) and synthetic spectra (continuous and dashed lines) of the early-R star HIP 84 266 around the A 
8015 A (top), and the A 4740 A C2 band regions (below), respectively. In both spectral ranges the synthetic spectrum giving the best 
fit (continuous line) was computed with C/N/O = 8.78/8.40/8.56 and 12 C/ 13 C= 7. The position of some 13 C 14 N lines used to derive 
the carbon isotopic ratio in the 8015 A region are indicated. Dashed lines in both panels correspond to synthetic spectra computed 
with 12 C/ 13 C = 4 and 12 C/ 13 C = 100, respectively. 
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~ ±0.3 dex in the gravity. Obviously, the most important source 
of uncertainty in the derived gravity is the luminosity. We esti- 
mate a maximum uncertainty in log g of ~ +0.9 dex due to the 
combined uncertainty in the mass, effective temperature and lu- 
minosity. However, the theoretical spectra are very sensitive to 
such a variation in gravity, so the actual uncertainty in gravity 
is actually lower than this maximum value (~ +0.5). The mean 
difference between the estimated gravity and the final adopted 
value in the analysis was ~ +0.1 dex for early-R stars. For late- 
R stars, the differences are below 0.5 dex except in the case of 
HIP 108205 (0.9 dex). 

We adopted initially a solar metallicity for a ll the stars acco- 
rding to the previous analysis by Dominy ( 1984). The r eference 
solar abundances are those from Aspl und et al.1 (120051) . The fi- 
nal metallicity of the star ([M/H]) was obtained as the average 
derived from a number of Fe, Ni, Mn and Zn lines. For solar- 
metallicity and moderate metal-poor stars these species scale 
approximately with Fe (i.e. [X/Fe]~ 0.0), thus, these elements 
can be safely used as metallicity indicators. The metallicity de- 
duced from the [X/H] ratios (with X — Fe, Ni, Mn and Zn) agreed 
within ±0.05 for the early-R stars and ±0. 10 for the late-R stars. 
A microturbulence parameter of 2 km s _1 was adopted in all the 
stars, which is a typical value for giant stars. This initial value 
was modified in the chemical analysis by fitting the profile of 
metallic lines apparently free of molecular blends. We note that 
the microturbulence velocity giving the best fit to the observed 
spectrum may change slightly with the wavelength range. The 
average microturbulence velocities are given in Table [3] The 
macroturbulence velocity adopted was between 5-10 km s _1 . 
Finally, the theoretical spectra were convolved with a gaussian 
distribution to mimic the intrumental profile. The final FWHM 
values adopted were within ~ 200 - 400 mA, depending on the 
wavelength range and the spectral resolution. 

In the analysis of carbon stars a critical issue is the determi- 
nation of the C, N, and O abundances since the actual C/O ratio 
in the atmosphere determines the dominant molecular bands and 
thus the global spectrum. The actual N abundance has a minor 
role in this sense. Unfortunately the O abundance is very diffi- 
cult to measure in carbon stars at visual wavelengths. Only in 
two early-R stars in our sample, HIP 39 118 and HIP 62 944, 
it was possible to derive the O abundance from the A 6300.3 
A [OI] line. For the other stars, the O ab undance was sc aled 
with the average metallicity. We note that iDominvl (I 19841) de- 
rived almost scaled to solar O abundances from molecular lines 
in the infrared in his sample of early-R stars. Nevertheless, the 
effect of the absolute O abundance on the theoretical spectrum 
is secondary with respect to the actual C/O ratio. This is because 
for a given C/O ratio (or difference [C/H] - [O/H]), there is a 
range of absolute O abundances (withi n a factor of three) giv - 
ing almost equal theoretical spectra (see lde Lavernv etaT1 l2006). 
To estimate the C and N abundances, we proceeded as follow: 
first, the carbon abundance was derived from the C2 Swan band 
lines at AA 4730 - 4750 A taking the nitrogen and oxygen abun- 
dances scaled with the average metallicity (except in the case of 
O for the two early-R stars mentioned above). The intensity of 
these C2 lines are not very sensitive to the N and O abundances 
adopted. From this region we also derive the 12 C/ 13 C ratio from 
some available 12 C 13 C and 13 C 13 C lines. Next, we synthesize the 
region at AA 8000 - 8050 A, which is dominated by the red sys- 
tem of CN molecule, to estimate the N abundance adopting the 
C abundance from the previous step. Then, we returned to the 
AA 4730 - 4750 A region to check the consistency of the theore- 
tical and observed spectrum with the new N abundance derived. 



Table 3. Atmospheric parameters and C, N, O abundances 



Star 


7-eff (K) 


l°g £ 


^(kms- 1 ) 


[M/H] 


C 


N 


O 


Late-R stars 


HIP 35 810 


3300 


0.0 


2.8 


-0.38 


8.40 


7.30" 


8.35 


HIP 36 623 


3300 


0.0 


2.1 


-0.27 


8.35 


7.70 


8.31 


HIP 91 929 


3300 


0.0 


2.5 


0.00 


8.68 


7.78" 


8.66 


HIP 108 205 


2500 


0.0 


2.3 


0.00 


8.67 


7.78" 


8.66 


HIP 109 158 


2500 


0.0 


2.5 


-0.02 


8.68 


7.78" 


8.66 


Early-R stars 


HIP 39 118 


4250 


2.0 


1.8 


-0.29 


8.77 


8.00 


<9.05* 


HIP 44 812 


3950 


1.5 


3.0 


-0.03 


8.75 


8.70 


8.66 


HIP 53 832 


4500 


2.5 


2.3 


-0.77 


8.08 


8.10 


8.06 


HIP 58 786 


4250 


2.0 


2.2 


-0.29 


8.63 


7.60 


8.47 


HIP 62 944 


4300 


2.4 


2.0 


0.12 


8.80 


8.30 


9.05* 


HIP 69 089 


4750 


1.5 


2.0 


-0.17 


8.40 


8.40 


8.46 


HIP 74 826 


4750 


2.0 


1.5 


-0.30 


8.43 


8.18 


8.36 


HIP 82 184 


4500 


2.0 


3.0 


-0.15 


8.42 


8.25 


8.51 


HIP 84 266 


4750 


2.0 


2.4 


-0.10 


8.78 


8.40 


8.56 


HIP 85 750 


3800 


1.2 


2.0 


-0.48 


8.38 


7.98 


8.36 


HIP 86 927 


4700 


2.4 


2.3 


-0.05 


8.47 


7.98 


8.66 


HIP 87 603 


4100 


2.0 


2.3 


-0.52 


8.60 


7.68 


8.36 


HIP 88 887 


3950 


1.5 


3.0 


-0.09 


8.80 


8.70 


8.66 


HIP 94 049 


4100 


2.0 


2.3 


-0.62 


8.48 


7.98 


8.36 


HIP 95 422 


3950 


2.0 


3.0 


-0.26 


8.73 


7.78" 


8.36 


HIP 98 223 


3800 


1.5 


1.9 


-0.79 


8.05 


7.10 


7.97 


HIP 113 150 


4500 


2.0 


2.4 


-0.47 


8.85 


7.60" 


8.36 



a N abundance scaled with the stellar metallicity. 

b O abundance derived from the A 6300.3 A [O I] line using the atomic line parameters 
from lCaffau etafl t20"08l) . 



This procedure was repeated several times until convergence was 
found. Typically two or three iterations were needed. For a given 
star, the maximum difference between the C and N abundances 
derived from the red and blue parts of the spectrum was less than 
0.2 or 0. 1 dex for the early- and late-R stars, respectively. The fi- 
nal carbon isotopic ratio obtained is the average value between 
those obtained from the fits to the A ~ 8000 A and A ~ 4750 A 
regions. The final C, N and O abundances together with the other 
atmospheric parameters are shown in Table [3] 

Synthetic spectra in LTE were computed fo r a specific star 
using the version 7.3 of the Turbospectrum code dAlvarez &Plez 
1998). We used the same extensive set of atomic and molecu- 
lar lines as in Papers II and III and refer to these works for de- 
tails. The list of lines is available from the authors upon request. 
Figures [3] to [5] show examples of theoretical fits in four differ- 
ent spectral regions used in the chemical analysis. The lithium 
abundance was derived using the resonant Li I at A 6707.8 A. 
This spectral range is also interesting because of the presence of 
one useful Y I line at A 6687.5 A. For most of the stars in the 
sample, the Li abundance has to be considered only as an upper 
limit. This is due to the severe blending in this spectral region, 
as can be seen clearly in Fig. [4] s-element abundances were de- 
termined mainly using the spectral window at AA 4750 - 4950 A 
(see Papers II and III) where several lines of Sr, Y, Zr, Ba, La, 
Nd and Sm are present (see Fig. [5]). For early-R stars we also 
used the Ba I-II lines at A 6498.9 and A 6496.9 A, respectively. 
We attempted to detect the radioactive element technetium from 
the resonant lines at A 4260 A and the recombination Tc I line at 
A 5924.5 A. We fail to detect the bluer line in all the stars mainly 
because of the low S/N ratio achieved in the spectra in that re- 
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Fig. 5. Observed and synthetic spectra in the A 4800 A region for the early-R star HIP 84 266 (top panel) and the late-R star HIP 
36 623 (bottom). Symbols as in Fig. [3] For each star two theoretical fits are shown: continuous line (best fit) is computed assuming 
[s/M] ~ 0.0, while dashed line with [s/M] ~ +0.5 in the star HIP 84266. For HIP 36 623, continuous line is computed with [s/M] 
~ +1.3 and dashed line with [s/M] ~ 0.0. 
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Species 


T e ff a 


logg 


[M/H] 




macro. 


C/O 


l2 C/ li C 


Total 


X 


+ 150/250K 


±0.5 


±0.2 


±0.5 km s" 1 


±1.5 km s _1 


±0.1 


±10 




Late-R stars 


Li 


+0.40 


±0.05 


±0.20 


+0.00 


±0.00 


+0.20 


±0.00 


±0.50 


Fe 


+0.30 


±0.00 




+0.05 


±0.10 


+0.10 


±0.00 


±0.35 


Rb 


±0.20 


±0.10 


±0.20 


+0.00 


±0.00 


+0.15 


+0.10 


±0.35 


Sr 


±0.15 


±0.00 


±0.10 


+0.05 


±0.20 


+0.10 


+0.05 


±0.30 


Y 


±0.15 


±0.00 


±0.10 


+0.20 


±0.25 


+0.05 


±0.00 


±0.35 


Zr 


±0.05 


±0.00 


±0.05 


+0.20 


±0.25 


+0.00 


±0.00 


±0.35 


Ba 


±0.15 


±0.10 


±0.10 


+0.05 


±0.00 


±0.00 


±0.00 


±0.25 


La 


±0.05 


±0.20 


±0.10 


+0.10 


±0.20 


+0.05 


±0.00 


±0.35 


Nd 


±0.05 


±0.20 


±0.05 


+0.20 


±0.20 


+0.05 


±0.00 


±0.35 


Sm 


±0.05 


±0.10 


±0.10 


+0.25 


±0.20 


+0.00 


±0.00 


±0.35 


Early-R stars 


Li 


±0.10 


±0.00 


±0.10 


+0.00 


±0.00 


+0.10 


±0.05 


±0.30 


Fe 


±0.05 


±0.00 




+0.05 


±0.15 


+0.10 


±0.00 


±0.20 


Rb 


±0.05 


±0.00 


±0.00 


+0.00 


±0.05 


+0.05 


±0.00 


±0.20 


Sr 


±0.10 


±0.00 


+0.05 


+0.10 


±0.15 


+0.05 


+0.00 


±0.20 


Y 


±0.10 


±0.00 


±0.25 


+0.00 


±0.10 


+0.15 


±0.00 


±0.40 


Zr 


±0.10 


±0.00 


±0.00 


+0.00 


±0.10 


+0.05 


±0.00 


±0.30 


Ba 


±0.10 


±0.10 


±0.20 


+0.00 


±0.10 


+0.05 


±0.00 


±0.30 


La 


±0.05 


±0.20 


±0.10 


+0.00 


±0.15 


+0.05 


+0.00 


±0.30 


Nd 


±0.15 


±0.20 


±0.05 


+0.05 


±0.20 


+0.00 


+0.00 


±0.40 


Sm 


±0.15 


±0.20 


+0.10 


+0.00 


±0.15 


+0.00 


+0.00 


±0.35 



" See text for details. 




lines at ^17780.6, 7802.5, and 7807.9 A not very much blended 
by molecular absorptions that can be used to derive the average 
metallicity. 



Table [4] shows the uncertainty in the abundances derived 
for the different chemical species due to the uncertainty in the 
model atmosphere parameters. To do this a canonical model at- 
mosphere with parameters r e ff/log g/[M/H]= 4750/2.0/0.0 and 
3300/0.0/0.0 was used for early- and late-R stars, respectively. 
Due to the very different sensitivity of the theoretical spectra to 
changes in the effective temperature, changes in the abundances 
due to this parameter were calculated assuming variations by 
±150 and ±250 K for early- and late-R stars, respectively. 



0.2 
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Fig. 4. Observed and synthetic spectra in the region of the Li line 
at A ~ 6707.8 A for the early-R star HIP 58786. Symbols as in 
Fig. [3] Continuous line is the best fit with log e(Li)= 1.0, and 
dashed lines are theoretical spectra computed with log e(Li)= 
0.0 and 1.5, respectively. 



gion. Upper limits on the Tc abundance were set in the stars 
HIP 39 118 and HIP 109 158 by using the recombination line. 
This result is compatible with the presence of some s-element 
enhancement (see Tabled in these two stars. However, because 
we can barely reproduce the blend at A 5924 A, these Tc up- 
per limits have to be taken with caution. Finally, Rb abundances 
were determined using the Rb I line at A 7800.2 A. In this region 
there are two Ni I lines at A 7788.9, 7797.6 A and three Fe I 



The final uncertainty in the abundance is found by consider- 
ing quadratically all the uncertainties including the uncertainty 
in the continuum position (< 5%) and, when more than 3 lines of 
a given element were used, the abundance dispersion among the 
lines of this element as an additional source of uncertainty. This 
is indicated in the last column of Table [4] From this table it can 
be appreciated that the error in the absolute abundances derived 
are in the range A[X/H] = ±(0.2 - 0.5) dex, being larger for 
late-R stars. When considering the abundances ratios ([X/Fe]) 
relative to iron, the error is reduced if the variation on a given 
atmospheric parameter affects in the same way the element of 
interest. The results of the chemical analysis are summarised in 
Tables [5] and [6] The last four columns in Table [6] show the av- 
erage enhancement of light ^-elements, ls= (Sr, Y, Zr), heavy 
s-elements, hs= (Ba, La, Nd, Sm), the s-element index [hs/ls], 
and the total average s-element enhancement. 
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Table 5. Abundances derived and relevant abundances ratios. 



Name 


C/O 


12 C/ IJ C 


Li 


Mn 


Fe 


Ni 


Zn 


[M/H] 


[C/M] 


[N/M] 


Late-R stars 


HIP 35 810 


1.12 


65 


0.20 




7.07 


5.93 


4.20 


-0.38 


0.39 




HIP 36 623 


1.10 


23 


-0.20 




7.18 


5.92 




-0.27 


0.23 


0.19 


HIP 62401 






















HIP 91 929 


1.05 


58 


-0.50 




7.45 


6.23 




0.00 


0.29 




HIP 108 205 


1.02 


90 


-1.50 




7.45 






0.00 


0.28 




HIP 109 158 


1.05 


85 


-1.00 




7.43 


5.93 




-0.02 


0.31 












Early- 


■R stars 












HIP 39118 


>0.50 


13: 


0.85 


5.10 


7.16 


5.93 


- 


-0.29 


0.67 


0.51 


HIP 44 812 


1.23 


5 


<1.00 


- 


7.42 


6.23 


4.60 


-0.03 


0.39 


0.95 


HIP 53 832 


1.05 


24 


<0.60 


4.69 


6.68 


5.43 


3.80 


-0.77 


0.46 


1.09 


HIP 58 786 


1.45 


70 


<1.00 


- 


7.16 


5.90 


4.15 


-0.29 


0.53 


0.11 


HIP 62944 


0.56 


22 


2.60 


5.49 


7.57 


6.28 


4.70 


0.12 


0.29 


0.40 


HIP 69 089 


0.87 


19 


1.80 


5.25 


7.28 


6.04 


4.50 


-0.17 


0.18 


0.79 


HIP 74 826 


1.17 


20 


<0.50 


5.10 


7.15 


5.95 


4.30 


-0.30 


0.34 


0.70 


HIP 82 184 


0.81 


10 


<0.40 


5.24 


7.30 


6.08 


4.45 


-0.15 


0.18 


0.62 


HIP 84266 


1.65 


7 


<1.05 


5.29 


7.35 


6.13 


4.50 


-0.10 


0.49 


0.72 


HIP 85 750 


1.05 


22 


<0.30 


4.89 


6.97 


6.97 


<3.85 


-0.48 


0.47 


0.68 


HIP 86 927 


0.65 


7 


<0.48 


5.35 


7.40 


7.40 


<4.55 


-0.05 


0.13 


0.25 


HIP 87 603 


1.74 


9 


<1.10 




6.93 


6.93 


4.20 


-0.52 


0.73 


0.42 


HIP 88 887 


1.38 


5 


<0.50 




7.36 


7.36 


4.30 


-0.09 


0.50 


1.01 


HIP 94049 


1.32 


9 


<0.40 




6.83 


6.83 




-0.62 


0.71 


0.82 


HIP 95 422 


2.34 


6 


<0.50 




7.19 


7.19 




-0.26 


0.60 




HIP 98 223 


1.20 


16 


<0.00 




6.66 


6.66 




-0.79 


0.45 


0.11 


HIP 113150 


3.09 


9 


<0.55 




6.98 


6.98 




-0.47 


0.93 





: indicates uncertain value. 



4. Results and discussion 

4. 1 . Comments on particular stars 

The detailed chemical composition reported in Tables [5] and [6] 
shows that there are several stars that do not fulfil one of the 
main characteristic that defines a R star, i.e., a carbon star with- 
out s-element enhancements. HIP 39 118 is one of these stars, 
showing a clear s-element enhancement, [s/M] — 1-0.54. In fact, 
our best fit to the Tc I line at A 5924 A in this star gives log 
e(Tc)< 1.20. This, together with the s-element enhancement 
found would be an indication of this star being on the AGB 
phase. However considering that the luminosity of this star is 
too low to be on the AGB phase (see Table |2), and the difficulty 
in reproducing the A 5924 A Tc blend, this Tc upper limit has to 
be considered with extreme caution and thus, it would be risky to 
conclude that HIP 39 1 18 is an AGB star. HIP 39 1 18 has other 
chemical peculiarities. It is not a carbon star since its C/O ratio 
(~ 0.50) is considerably lower than unity even considering the 
uncertainty in the derivation of this ratio. Indeed, the A 6300 A 
[OI] line is very strong in this star which indicates some O en- 
hancement, [0/Fe]< 0.68. Titanium, another alpha element, is 
also enhanced in this star [Ti/Fe]= +0.3. Due to its low lumi- 
nosity (Mbd = 0.47) and ^-element enhancement one is tempted 
to ascribe this star as a classical barium star that are known to 
be all binaries. In the case of HIP 39 1 18, there are only two ra- 
dial velocity measu rements within a 5 days interval separation 
(Plata is et al.l l2003). With such short time span compared with 
the long orbital periods of the classical barium stars, no conclu- 
sion may be drawn regarding binarity. Moreover, the moderate 
Li abundance in this star, log e (Li)= 0.85, is difficult to explain 



in a mass transfer scenario. Thus, the spectral classification of 
this star is uncertain and we just put it in the group of chemi- 
cally anomalous red giant stars, but certainly not an early-R star. 

HIP 53 832 has an abundance pattern typical of a classic CH- 
type star: a moderate metal-poor carbon-rich ([M/H] = -0.77, 
C/O = 1.05) star with a relatively low carbon isotopic ratio, 
( 12 C/ 13 C= 24) and an important ^-element enhancement ([s/M] = 
+ 1.26). In fact, the binary nature of this star has been con- 
firmed by radial velocity measurements by Plat ais et"ail ([2003) 
and thus, there is no doubt about its classification as a CH star. 

HIP 62944 is one of th e first Li-rich giants discovered 
dWallerstein & Snedenl |1982|) . Our Li abundance in this star, 
log e(Li)= 2.60, agrees with that deduced by these authors. 
The CNO content and the metallicity derived here are also in 
agreement. Again this star is not a carbon star as the C/O ra- 
tio is near the solar value (see Table [5). The low isotopic ratio 
( 12 C/ 13 C= 22) derived here is similar to those found in others Li- 
rich giants (Ide La R eza 2006). The origin of Li in RGB stars is 
not well understood. Some kind of n on-standard mixing mecha- 
nism during the RGB p hase (see e.g. lPalacios et al.ll2001l 12006: 
Guandalini et al. 2009) seems to be needed. In any case, HIP 
62 944 is a K giant, not a R star. 

HIP 69089 presents some C enrichment although its C/O 
ratio is below unity (0.87, see Table |5). As in the case of HIP 
62 944, this star is also Li-rich (log e(Li)= 1 .80) with a low car- 
bon isotopic ratio (19 ) typical in these stars. Th e Li abundance 
previously derived by iLuck & Challenerl d 1995b in this star, log 
e(Li)= 2.04, is compatible with the value derived in the present 
work. The ^-element abundances are solar scaled. On the basis 
of these values we classify it as another Li-rich K giant. 
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Name 


[Rb/M] 


[Sr/M] 


[Y/M] 


[Zr/M] 


[Ba/M] 


[La/M] 


[Nd/M] 


[Sm/M] 


[ls/M] 


[hs/M] 


[hs/ls] 


[s/M] 


Late-R stars 


HIP 35 810 




0.26 


0.27 


0.29 


0.78 


0.55 


0.83 


0.17 


0.27 


0.58 


0.31 


0.43 


HIP 36 623 


0.10 


0.95 


1.31 


1.13 


1.72 


1.64 


1.72 


0.56 


1.13 


1.41 


0.28 


1.27 


HIP 91 929 






-0.01 


0.01 


-0.02 








0.00 


-0.02 


-0.02 


-0.01 


HIP 108 205 






0.59 




0.73 








0.59 


0.73 


0.14 


0.66 


HIP 109 158 






0.81 




0.85 








0.81 


0.85 


0.04 


0.83 


Early-R stars 


HIP 39118 


0.30 


0.22 


0.68 


0.70 


0.49 


0.46 


0.64 


0.58 


0.53 


0.54 


0.01 


0.54 


HIP 44 812 




-0.07 


0.03 


-0.26 


-0.04 


0.13 




0.12 


-0.10 


0.07 


0.17 


-0.01 


HIP 53 832 


<0.70 


0.85 


1.06 


0.78 


1.65 


<2.04 


1.42 


1.36 


0.90 


1.62 


0.72 


1.26 


HIP 58 786 




0.09 


-0.06 


-0.30 


0.02 


-0.01 




-0.01 


-0.09 


0.00 


0.09 


-0.04 


HIP 62 944 


-0.09 


-0.02 


-0.17 


-0.41 


-0.29 


-0.02 


-0.12 


0.07 


-0.20 


0.09 


0.11 


-0.15 


HIP 69089 


-0.13 


0.05 


0.16 


-0.12 


0.00 


0.07 


0.07 


0.06 


0.03 


0.05 


0.02 


0.04 


HIP 74826 


0.01 


-0.02 


0.02 


0.01 


0.08 


-0.04 


0.00 


-0.01 


0.00 


0.01 


0.00 


0.01 


HIP 82184 


-0.01 


0.00 


-0.07 


-0.14 


0.08 


0.15 


-0.10 


0.00 


-0.07 


0.03 


0.10 


-0.02 


HIP 84266 


-0.04 


0.00 


-0.21 


-0.11 


-0.07 


0.22 


0.05 


0.14 


-0.11 


0.09 


0.19 


-0.01 


HIP 85 750 


-0.05 


0.76 


0.97 


0.59 


1.56 


1.65: 


1.53 


1.27 


0.77 


1.45 


0.68 


1.11 


HIP 86 927 


0.18 


0.00 


-0.01 


0.31 


0.20 


0.32 


0.40 


0.19 


0.10 


0.28 


0.18 


0.19 


HIP 87 603 


0.05 


0.12 


0.11 


-0.07 


0.45 


0.52 


0.12 


0.11 


0.05 


0.30 


0.25 


0.18 


HIP 88 887 


-0.07 


-0.13 


-0.02 


-0.20 


-0.18 


-0.11 


-0.11 


-0.12 


-0.12 


-0.13 


-0.01 


-0.12 


HIP 94049 






0.02 




-0.05 








0.02 


-0.05 


-0.07 


-0.01 


HIP 95422 






-0.04 




0.11 








-0.04 


0.11 


0.15 


0.04 


HIP 98 223 




0.67 


1.13 


1.00 


1.25 




1.44 


1.18 


0.93 


1.29 


0.36 


1.11 


HIP 113 150 






-0.03 


-0.02 


0.00 






-0.04 


-0.02 


-0.02 


0.00 


-0.02 



Is = <Sr, Y, Zr), hs = <Ba, La, Nd, Sm> and s= <ls + hs >. 



We derive significant s-element overabundances in the car- 
bon star HIP 85 750 (see Table [6}. Radial velocity variations 
have been detected in this star (Table[TJ. Since it is a moderately 
metal-poor star, we classify it as a CH-type star. 

According to our analysis HIP 86927 is not a carbon star 
(C/0= 0.65), and shows no s-element enhancement. Its carbon 
isotopi c ratio is l ow (7 ) and agrees with the previous determina- 
tion by Dominvl (ll984l) . Unfortunately, this author did not derive 
any other abundance signature in this star. Further, the luminos- 
ity of this star (Mboi = 0.22) is too low for a typical early-R star 
(see Sect. 12.6) , thus we believe that its characteristics are closer 
to a normal K giant which has undergone some extra-mixing re- 
vealed by its low 12 C/ 13 C ratio. 

Finally, HIP 98 223 is a metal-poor ([M/H] = -0.79) carbon 
star showing a large enrichment in s-elements ([s/M]= +1.11). 
These are the typical abundance signatures of the CH stars. 

Among the late-R st ars we have two peculiar stars. HIP 
36623 is a symbiotic star (|Johnson et al.lll988l : lBelczv nski et al. 
2000: iMunari & Zwitterll2002l) whose luminosity (see Table 
and chemical composition (see Tables [5] and |6]) are compatible 
with those expected at the AGB phase. Perhaps the envelope of 
this star has undergone two periods of mixing: the first one by 
the material accreted from the primary star, and the second one, 
a self-contamination triggered by third dredge-up episodes on 
the AGB phase. Thus, it is probably a normal N-type star. The 
spectrum of this star looks similar to the other late-R in the sam- 
ple without any evidence of contamination by emission lines or 
emission continua coming from the hot UV radiation of the com- 
panion (probably a white dwarf). The broad emission features 
around 6825 and 7082 A (due to O VI photons Raman scattered 



off H I) are not present either. These features are found in some 
symbiotic stars (e.g. AG Dra, ISmith et id] [1996"). Thus, we be- 
lieve that the abundances of the s-elements in this star (mostly 
derived from the blue part of the spectrum) are not underesti- 
mated because of a contribution to the spectral continuum by 
the companion star. On the other hand, we do not find ^-element 
overabundances in the late-R HIP 91 929. However its spectrum 
is identical to those of N-type stars. In Paper II we found a few 
N-type stars where the s-element content was compatible with 
no enhancements within the error bars. Perhaps HIP 91 929 is 
another example of this. Maybe these stars are at the beginning 
of the thermally pulsing AGB phase (as their luminosity might 
indicate) and did not have time yet to pollute the envelope with 
i-elements. 

In summary, out of seventeen stars classified as early-R stars 
we have found seven stars (~ 40 %) that are not R stars. Since the 
stars in our sample were randomly chosen from the Hipparcos 
catalogue (the only condition was to be observable from the 
north hemisphere) and although the sample is limited, we con- 
firm the previous claims that among the R stars there is a sig- 
nificant number of stars which are wrongly classified. This is an 
important result since it might discard the previous belief (see 
Sect. Q]) that the R stars (namely early-R) represent a frequent 
stage in the evolution of low-mass stars. 

4.2. CIO and 12 CI 13 C ratios 

Fig. [6] shows the C/O and 12 C/ 13 C ratios found in our stars 
(C= 12 C+ 13 C). The stars reclassified here (see previous section) 
are indicated by open symbols (squares, CH stars and trian- 
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gles, K giants) In the figure are indicated the expected values 
of the 12 C/ 13 C ratio according to standard low-mass stellar mod- 
els (e.g. ICristallo et all 12009) after the first dredge-up, and the 
AGB phase, respectively. The expected range in the AGB phase 
has been computed considering the e ffect of an extra-mixing 
mech anism after the first dredge-up ( Bo othrovd & Sackmannl 
1999), which would reduce further the 12 C/ 13 C ratio to a typi- 
cal value of ~ 12. Such 1 2 C/ 13 C ratios are frequently found in 
Population I giant stars dGilrovlll989t iGilrov & Brown] 1 19911) 
and can be only explained ass uming a non standard mixing 
process (Char bonnel et al.lll998l) . It is evident from Fig. [6] that 
most of the true early-R stars have carbon isotopic ratios far 
below the minimum values expected after the first dredge-up 
(RGB phase) and the thermally pulsing AGB phase. Further, 
since early-R stars show also N enhancements (see Table [5] and 
lDominy|[l98ll) . the low 12 C/ 13 C ratios together with the N over- 
abundances (<[N/M]>= +0.6) would indicate that the material 
in the envelope of these stars has been processed through the 
CNO bi-cycle (almost at the equilibrium). If that were the case, 
material exposed to the CNO bi-cycle is expected to be depleted 
in 16 by a factor of at least ~ 40, the exact value depending on 
the burning temperature. Even considering the large uncertainty 
in the O abundance in our stars (actually we scale the O abun- 
dance with the average metallicity, see Sect. 3), we can safely 
discard a large O depletion from the analysis: the C/O ratios 
would be much higher than the values found here and because 
few C would be locked into CO, the C-bearing molecules (CN, 
CH and C2) would appear very strong in the spectrum which is 
not observed in any of the early-R stars. Furthermore, fo r the few 
early- R stars where O abundance has been determined (Dominv 
1984) no evidence of oxygen depletion is found. This fact casts 
some doubt on the idea that the mixing of CNO-equilibrium ma- 
terial is the only responsible for the early-R stars. The C and N 
enhancements found are rather evidence that the material that we 
currently observe in the envelopes of early-R stars has been pro- 
cessed by both the CN cycle and He-burning. Probably, C-rich 
material was mixed with protons on a time scale short enough 
and at a temperature low enough (T< 7 x 10 7 K) to allow only 
the reactions 12 C(p,y) I3 C and 13 C(p, y) 14 N, without any signif- 
icant operation of the ON-cycle. In any case, the determination 
of the oxygen isotopic ratios in early-R stars would be a valu- 
able tool to evaluate the role played by the CNO bi-cycle in the 
surface composition of these stars. 

On the other hand, the I2 C/ I3 C ratios derived in late-R stars 
are significantly higher than those for the early ones. In these 
stars observations in general agree with AGB model predic- 
tions (see Fig. |6). An exception is the star HIP 36623 with a 
12 C/ 13 C ratio below the predicted values for the AGB carbon 
stars. However, thes e low 12 C/ 13 C values have been found a lso in 
many N-type stars dOhnaka & Tsui ill 19961: [Abia et al.ll2002l) and 
are currently interpreted as an evidence of the operation of a non 
standard mixing process also in low-mass stars during the AGB 
phase. Differences between the early- and late- R stars are also 
evident in the C/O ratio. Early-R stars show a wider range in the 
C/O ratio while late-R stars are concentrated ver y close to 1 , as 
in the galactic N-type s tars of similar metallicity ( Lambert et al.l 
1986: [Abia et al.ll2002l) . This is consistent with late-R stars hav- 
ing higher mass than early-R stars because it is more difficult to 
increase the C/O ratio of a more massive star. 

4.3. Lithium 

As noted in Sect. [3] most of the Li abundances derived in early- 
R stars should be considered upper limits. In cooler late-R stars, 
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Fig. 6. C/O vs. 12 C/ 13 C ratios. Symbos as Fig. 2. The arrows in- 
dicate the expected range in the 12 C/ 13 C ratio according to the 
standard evolutionary models of low-mass stars in the RGB and 
AGB phases (see text for details). 
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synthetic spectra are much more sensitive to changes in the Li 
abundance (but note however the large formal uncertainty in the 
Li abundance, see Table 4). Fig. Q presents the Li abundances 
derived in late- and early-R stars (excluding the stars reclassi- 
fied in the previous section) and the luminosity. We remind that 
the typical stellar mass of early- and late-R stars is probably very 
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Fig. 8. Comparisons of the derived s-element ratios ([X/Fe]) in 
the late-R star s HIP 36 623 and HIP 35 810 with the theoretical 
predictions bv lCristallo et al.l (120091) for two different low-mass 
AGB models. In each panel is indicated the stellar mass, metal- 
licity and specific thermal pulse that provide the best fit to the 
observed abundances. 



different (see Sect. 12. 31 thus, one should not establish a sequence 
of Li depletion in R stars from this figure. Instead, the appar- 
ent Li abundance vs. luminosity relationship have to be inter- 
preted separately. By comparing stars of the same spectral type 
we found that early-R stars present very similar Li abundances 
(in the range log e(Li)~ 0.5 - 1 .0) without any evident correlation 
with the luminosity. These upper limits in the Li abundances are 
higher than the expected values in post-tip RGB stars (marked 
with the dotted line at log e(Li)= 0.0 in Fig. 7, e.g. ICastilho 
2000). Although, it would be necessary to confirm the possible 
Li detection with higher resolution spectra, this is of some rele- 
vance because it may constrain the scenarios proposed to form a 
R star. 

In late-R stars, at contrary, Li abundances decrease with the 
increasing luminosity as one would expect to occur as the stellar 
envelope goes deeply toward the interior of the star during the 
ascend along the AGB phase. Indeed, the Li abundances derived 
in our late- R stars are similar to those derived in normal AGB 
carbon stars dDenn et al ,1ll99lHAbia et aljl993l) . which supports 
again the idea that they are identical to the N-type stars. 

4.4. Abundances of heavy elements 

Table [6] summarises the heavy element abundances derived in 
our stars. In the sample of early-R stars, and excluding those 
stars that we have reclassified as belonging to other spectral 
types (see above), it is evident that these stars do not have s- 
element overabundances [X/M]~ 0.0 within the error bars. This 
is the case of HIP 44 812, HIP 58 786, HIP 74 826, HIP 82 184, 
HIP 84266, HIP 87 603, HIP 88 887, HIP 94049, HI P 95422 
and HIP 113150. Thus, we confirm the analysis by |Dominv 



(1 19841) in this kind of stars. Late-R stars do have s-element 
enhancement^ at the same level than those found in galactic 
N-typ e stars of similar metallicity, <[s/M]>= 0.64 (lAbia et al.l 
2002). In fact, we detect Tc (using the weak recombination line 
at A 5924.47 A) in HIP 109 158 and it may be also present in 
HIP 108 205. This is usually interpreted as evidence of the in 
situ production of the s-elements in the star, the envelope being 
polluted by the operation of the third dredge-up during the AGB 
phase. These stars are named intrinsic. Fig.|8]shows an example 
of a detailed reproduction of the observed heavy element pattern 
in two late-R stars by theoretical s-process nucle osynthesis cal- 
culati ons in low-mass AGB stars, (see details in ICristallo et al.l 
2009). Each panel indicates the stellar mass model and thermal 
pulse giving the best fit to the observed pattern. The metallicity 
in the model (Z) was chosen as close as possible to the observed 
metallicity in the stars. By selecting a stellar mass model, metal- 
licity and specific thermal pulse, fits of similar quality can be 
found for the other late-R stars in our sample, as well as for the 
early-R stars reclassified as CH-type (HIP 98223, HIP 85 750 
and HIP 53 832) which also show heavy element enhancements. 
In this case, we can mimic their extrinsic nature (i.e. the chem- 
ical composition of the envelope is due to the pollution by the 
matter accreted from a companion) including an extra-parameter 
in the model: the dilution factor, defined as the ratio between the 
mass of t he envelope and the mass transferred from the primary 
star (e.g. lBisterzo et al1 l2006). The dilution factor needed are in 
the range 2-3. These values are consistent with a secondary 
star having a large envelope, most of the CH-type star are indeed 
giants. 

We close this section realising once more that the most pro- 
bable neutron source producing the observed s-element pat- 
tern in these stars is the 13 C(a, n) lfi O reaction which is found 
the main neutron donor in low-mass AGB stars (Lambert et al. 
1995; I Abia et all 12001 1) . This is concluded from the observed 
[Rb/Sr, Y,Zr] < (see Table|6]l. In AGB stars of intermediate mass 
(M > 3 M©) probably dominates the 22 Ne(a, n) 25 Mg neutron 
source resulting in higher n eutron densities (N„ > 10 10 c m" 3 ) 
and so [Rb/Sr, Y,Zr]> (e . g. iGarcia-Hernandez etaT1l2006l) . 



4.5. Evolutionary status 

We have shown in the previous analysis that the true early- and 
late-R stars have different properties. Most of the late-R stars 
are long -period variables of the SR or Mira types (see Table [TJ, 
and they often exhibit excess emission at 12 fim due to dust, 
indicative of mass loss. In the colour-colour diagram they oc- 
cupy the same location as N-type stars and have similar lumi- 
nosities (or slightly lower). We have also shown that there is no 
significant difference in the chemical composition of late-R stars 
and N-type stars, both types showing s-element enhancements 
at the same level. They are thus closely related, or are identical. 
Definitely late-R stars are AGB stars and, most probably, mark 
the bottom of the AGB phase. 

Concerning early-R stars, our study (although limited in the 
number of stars) confirms previous claims that most carbon stars 
catalogues contain a significant number of stars wrongly classi- 
fied as early-R stars. In a sample of seventeen early-R stars we 
have found seven (~ 40%) objects that are K giants or CH-type 
stars. This finding, together with the definite ascription of late- 
R stars as N stars certainly reduce the number of real (early) R 



4 The late-R star HIP 62401 is excluded from the chemical analysis 
because its spectrum shows very broad lines that we cannot reproduce. 
This is typically found in Mira variables, see Table[T] 
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type stars among all the giant carbon stars. Therefore, R stars 
are indeed scarce objects, but we still lack an evolutionary sce- 
nario to explain their observational properties. We have to con- 
sider at least three possibilities: i) mass transfer in a binary sys- 
tem, ii) original pollution, and iii) a non-standard mixing mech- 
anism able to mix carbon to the surface. Due to their position in 
the H-R diagram, and their peculia r chemical composition, the 
favou red scenario is the later one dDominvlll984t iKn app et al. 
20011) triggered by an anomalous H e-flash (e.g. lMengel & Gross 
197fj|Paczvnski & Tremaindl 19771) . 

As noted in Sect. 1, all early-R stars seem to be single stars: 
so far, attem ps to detect radi al velocity variations due to binarity 
have failed (lMcClure1ll997h . This is indeed a very improbable 
figure since in any stellar population one would expect a mini- 
mum of ~ 30% binary systems. McClure (1997) used this argu- 
ment to suggest that R stars were initially all binaries and that 
they coalesced into a single object during their evolution. In this 
hypothesis one might argue that the carbon enrichment that we 
currently observe was a consequence of mass transfer prior to 
the coalescence. However, it is extremely difficult to form car- 
bon s t ars from mass transfer at near solar metallicity (Abi a et all 
120031: iMasseron et alj|2009l alt hough a few exce ptions seem to 
exist, e.g., BD +57° 2161, see IZa£s etai1l2005l) . thus there is 
very little or not space for the mass transfer scenario. The se- 
cond scenario (carbon excess in the original gas cloud) it also 
ruled out because one should expect to find carbon stars of near 
solar metallicity at earlier evolutionary stages (main sequence, 
turn-off, sub-giants stars etc.). This is the case of the carbon rich 
metal-poor stars discovered in the galactic halo (e.g. Aok i et ail 
120071: IMasseron et alj[2009b : they are found along all the H-R 
diagram from the main sequence to the RGB. However, no a 
sole star with similar characteristics to those of early-R has been 
found in earlier evolutionary stages. In the same line, we might 
ask where are the descendants or early-R stars? Once leaving 
the He-core burning phase, they will evolve and ascend the giant 
branch for a second time, i.e. they should become AGB stars. 
The carbon stars of J-type have be en traditionally pr oposed as 
the daughters of early-R stars (e.g. lLlovd Evanslll986l) because 
of their chemical similarities (both are carbon stars without s- 
element enhancements) and the higher luminosity of the J-stars 
(typical of the AGB phase). The large Li abundances typically 
found in J-type stars (~ 80% show log e(Li)> 0.5 - 1.0, Paper 
I) are consistent with the Li excess that we have found here in 
early-R stars, supporting this idea. Nevertheless, this picture is 
far from clear since, for instance, the galactic distribution of the 
J-type stars is different from that of the early-R stars: J-type 
stars are located mainly in t he galactic thin d isk (only ~ 17% 
of the J-type stars studied by Chen et aT1l2007l have galactic lat- 
itude values |b| > 25°, typical of the thick disk stars). Secondly, 
a few J-type stars have been found in bina ry systems (e.g. BM 
Gem, EU And, UV A ur and UKS-Cel, iBarnbaumet all 1 199 It 
Belczvnski et al. 2000). Further work is necessary to determine 
if early-R stars are, in fact, the progenitors of some of the J-type 
stars. 

The solution seems to be related with the third sce- 
nario during or after the He-flash. However, as commented in 
Sect. [T| one-dimensional (ID) attempts and very recent three- 
dimensional (3D) numerical simulations of the He-flash for 
stars with a metallicity close to solar, as observ e d in early-R 
stars ([Dearborn etal.ll2006l lLattanzio etaLI l2006t iMocak et all 
2008albir have fail to mix carbon with the envelope. Binary 
mergers have been suggested as the mechanism able to provoke 
an anomalous He-flash. In the merger, a fast rotating He core 
is supposed to be formed with, as a consequence, a strong off- 



center He-flash. This stronger He - flash e ventually would mix 
carbon to the surface. Izzard et al] d2007l) explored statistically 
the different binary scenarios that might produce a R star. They 
discussed actually eleven possible channels, the most favourable 
cases being the merger of a He white dwarf with a RGB star 
and that of a RGB star with a Hertzsprung gap star. In fact, 
their binary population study produces ten times as many stars 
as required to match the early-R to red clump ratio observed. 
This discrepancy (that becomes more severe considering our 
results) was, however, interpreted by these authors as a posi- 
tive result since they expect that only a very small fraction of 
the favourable systems would ignite He while rotating rapidly 
enough to provoke the conditions for carbon mixing into their 
envelope. We have tested this scenario by parametrised ID sim- 
ulations of the merging and, for the first dynamical part, by 3D 
Smoothed Particle Hydrodynamics (SPH) simulations. An ex- 
tended discussion of these numer ical simulations wi ll be pre- 
sented in an accompanying paper (Piersanti et al. 2009), but we 
can advance here that in the most favourable merging channels 
according to Izzar d et al. (2007), we do not obtain any carbon 
mixing (see also lZamorall2009l). 

Very recently, Wallerstei n et alj ((2009) have found several 
carbon-rich (C/O > 1) RR Lyrae stars with similar N enhance- 
ments to those in early-R stars and have suggested that carbon 
might have mixed to the surface during the He-flash by an still 
unknown physical conditions. More theoretical work on this sub- 
ject is needed to determine such conditions. 

5. Summary 

We have shown that the early- and late-type R stars have dif- 
ferent properties including their chemical composition. Late-R 
stars have almost identical chemical figures than normal (N- 
type) AGB carbon stars, they occupy a similar position in the HR 
diagram, are long period variables and present infrared excesses 
due to dust. These stars probably mark the bottom of the AGB 
phase and can be considered thus identical to the N stars. For the 
early-type R star s we confirm th e chemical features found in the 
early analysis by Dominyj d 19841) and the fact that many CH-type 
and/or K giants have been erroneously classified as early-type R 
stars in the available carbon stars catalogues. So early-type R 
stars are indeed rare objects accounting a much smaller fraction 
among the giant carbon stars than previously thought. The lo- 
cation of these stars in the red clump provides strong support 
for their formation through an anomalous He-flash that would 
mix carbon to the surface. However, all the theoretical attempts 
made up to date have fail to reproduce such mixing. The origin 
of early-type R stars still remains a mystery. 
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